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The "pushing gate" proposed by Cirac and ZoUer for quantum logic in ion traps is discussed, 
in which a force is used to give a controlled push to a pair of trapped ions and thus realize a phase 
gate. The original proposal had a weakness in that it involved a hidden extreme sensitivity to 
the size of the force. Also, the physical origin of this force was not fully addressed. Here, we discuss 
the sensitivity and present a way to avoid it by choosing the spatial form of the pushing force 
in an optimal way. We also analyse the effect of imperfections in a pair of tt pulses which are 
used to implement a "spin-echo" to cancel correlated errors. We present a physical model for 
the force, namely the dipole force, and discuss the impact of unwanted photon scattering, and of 
finite temperature of the ions. The main effect of the temperature is to blur the phase of the gate 
owing to the ions exploring a range of values of the force. When the distance scale of the force 
profile is smaller than the ion separation this effect is more important than the high-order terms in 
the Coulomb repulsion which were originally discussed. Overall, we find that whereas the "pushing 
gate" is not as resistant to imperfection as was supposed, it remains a significant candidate for ion 
trap quantum computing since it does not require ground state cooling, and in some cases it does not 
require the Lamb-Dicke limit, while the gate rate is fast, close to (rather than small compared to) 
the trap vibrational frequency. 

PACS numbers: 03.67.-a, 42.50.-p 



I. INTRODUCTION 



Since the first recognition of the advantages of ion 
traps for quantum computing [1], various proposals to 
achieve quantum gates between pairs of trapped ions 
have been put forward [2] . In all cases each ion stores one 
or more qubits in its internal state, and most proposals 
have envisaged two or more ions in the same harmonic 
well, with the joint motional degree of freedom (normal 
mode of oscillation) serving as a further carrier of quan- 
tum information, which can be coupled to any chosen 
ion by laser excitation [3-6]. More recently a method 
of a qualitatively different form was proposed, in which 
the ions need not be in the same harmonic well, and 
a quantum gate is achieved through a more direct use 
of the Coulomb repulsion between ions [7]. Two ions 
in neighbouring harmonic wells are pushed by a force 
which depends on their internal state. Therefore, their 
mean separation depends on their internal state while 
the force acts. It is shown in Ref. [8] that the phase ac- 
quired from the Coulomb energy has the form of single- 
qubit rotations combined with the controUed-phase gate 
|00) (00|-h|01) (01|-|-|10) (lOl-fe*'' 111) (11|. This "push- 
ing gate" is significant for three reasons. First, it was 
shown [8] that in the right conditions the gate is unusu- 
ally insensitive to the details of the motion of the ions, 
i.e. the quantum phases depend primarily on the mean 
positions of the ions, and are insensitive, for example, to 
the ion temperature. Secondly, the gate can be faster 
than those previously proposed. Thirdly, the fact that 
the gate can act on ions in separate wells is useful for 



scaling the computer up to large numbers of ions, be- 
cause it is possible to move the ions around without ever 
having to "tease apart" two ions which were initially in 
the same well. 

For these reasons the "pushing gate" is a promising 
idea, but various aspects were left unclear in work up till 
now, as follows. 

(i) The method generates single-qubit rotations simul- 
taneously with the desired two-qubit phase gate. In 
the original work [8] it was assumed that these rota- 
tions would be unproblematic, since they are known 
from the experimental parameters and can be undone 
afterwards, if so desired, by any available single-qubit 
gate method. Therefore, the explicit calculation of these 
phases was not carried out. However, if the rotation an- 
gles are large compared to tt they will be sensitive to 
experimental imprecision. Therefore, it is necessary to 
calculate them and if possible, find ways to suppress or 
compensate for them. 

(ii) A state-dependent pushing force was introduced 
without analysing in sufficient detail the physical mech- 
anism which gives rise to the force. 

(iii) The proposal assumed that the force is indepen- 
dent of position, and it was unclear to what extent 
the good behaviour of the gate relied on that assump- 
tion. 

(iv) The effect of imprecision or fluctuations in the pa- 
rameters was not considered. 

The reason it is important to clarify these points is 
that until this is done, it is not clear whether the good 
performance of the gate is genuine, or whether it is an ar- 
tifact of some unphysical assumption, such as an unrea- 
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sonably high laser intensity or an unreasonably high pre- 
cision in parameters. This paper will address the points 
just listed. 

In Sec. II we introduce the basic concept of a gener- 
alised phase gate in order to bring out the physics without 
reference to any specific physical system. We also exam- 
ine the double 7r-pulse "spin-echo" which was suggested 
in Ref. [8] as a way to improve the fidelity. The "spin- 
echo" is used to cancel a correlated error which might 
be large, so there is a danger that a small imprecision in 
this part of the evolution might result in a large loss of 
fidelity. Therefore, we examine the influence of impreci- 
sion in the n pulses. 

In Sec. Ill we introduce the ion trap system. First, 
we show how the gate works in a qualitative way, 
and then we write down the Hamiltonian and deduce 
the classical trajectories of the ions under the influence 
of the two trapping potentials, their Coulomb repulsion, 
and the state-selective and time-dependent pushing force 
which is, for this section, assumed to be independent of 
position. 

Sec. IV then obtains all the quantum phases which 
arise in the evolution in the semiclassical approxima- 
tion, where the trajectories are classical and the quantum 
phase is obtained from a path integral along the classical 
trajectory. Our main purpose is to obtain the single- 
qubit rotation phases which were not calculated in 
the original paper [8]. In order to make the treatment 
self-contained we also rederive the two-qubit phase which 
is the central feature of the gate. 

Sec. V considers an explicit choice for the pushing 
force. The most natural choice is the dipole force from 
a non-resonant laser beam. This has the unavoidable con- 
sequence that unwanted scattering of photons will occur, 
which reduces the fidelity. We calculate the photon scat- 
tering and hence the loss of fidelity. 

In Sec. VI we examine the influence of fluctuations in 
the laser intensity. This is important because the single- 
qubit phases introduced by the gate are large, so that 
in some parameter regimes small changes in laser inten- 
sity can cause single-qubit rotations much larger than tt, 
destroying the fldelity. We propose a way to largely cir- 
cumvent this problem by appropriately choosing the di- 
rections of the forces and the distance scale of the poten- 
tial energy function (e.g. AC Stark light shifts), which 
provides the force. 

In Sec. VII we examine the effect of position depen- 
dence of the force, which is a major constraint on the tem- 
perature of the ions. It is also non- negligible since tightly 
focused laser beams, or a laser standing wave, must be 
used to get sufficient dipole force without unreasonably 
high laser power. 

In Sec. VIII we consider the total fidelity of the gate, 
including the effects of imperfect tt pulses, photon scat- 
tering, positional dependence of the force, and non-zero 
temperature of the ions. We give example values for 
the calcium ion, and discuss the results in comparison 
to other gate methods in ion traps. 



II. GENERAL CONCEPT 

In this section we shall introduce a general concept 
of a two-qubit phase gate. We discuss the fidelity of 
the phase gate and also the fidelity of the gate incorpo- 
rating a double 7r-pulse or "spin-echo" method to cancel 
one source of infidelity. Finally, we treat the latter case 
with imperfect tt pulses. 

A. Phcise gate 

Consider a system of two interacting qubits described 
by the evolution operator G in the computational basis 
{|00), |01), 1 10), 1 11)} having the form 

1 

a, 0=0 

= diag{e'^""^°°,e*®°\e*®i°,e*®"}, (1) 

where the phases Qaf3 have different values in general 
because the states \a(3) = \a)i (8) |/3)2 can have a differ- 
ent energy during the interaction and hence a different 
action integral. A physical system which satisfies the evo- 
lution (1) will be introduced and described in Sec. III. 

The operation G can be transformed into a two-qubit 
phase gate which generates a phase i? if and only if 
both qubits are in the logical state |1), that is 

= diag{l,l,l,e'''}. (2) 

The transformation G ^ P^ can be accomplished by 
local operations (single-qubit rotations) applied on each 
qubit 

S = Si®S2, (3) 

with 

5i = ^1(610-600), (4a) 
^2 = ^2(601-600), (4b) 

where 

Z,-(e)^|0),(0|e'®/2 + |l),-(l|e-'®/2 (5) 

is a rotation by the angle O about the z axis of the Bloch 
sphere. Bringing together Eq. (1) and (3) we can realise 
the phase gate (2) in the following way 



100) 


e"^^"" |GG) ^ 


|00), 


101) 


— > e'^°^ |01) — > 


101), 


|10) 


— > e*®i° |10) — ^ 


|io), 


111) 


— > e'®" |11) — > 


e'''|ll). 



(6) 



where 

= en - eio - eoi + Goo , (7) 
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and a global phase (8oi+6io)/2 was neglected. The two- 
qubit phase gate with the phase i? = tt is equiva- 
lent (up to additional single-qubit rotations) to a two- 
qubit controlled- NOT (CNOT) gate, which together with 
single-qubit rotations forms a universal set of gates for 
quantum computation [9]. 



B. Imprecision of the phcise gate 

In practice for a particular physical system the phases 
Qaf3 in Eq- (1) arc not accurately known and they may 
vary from one realisation of the gate to another. On 
the other hand, we can assume the single-qubit rota- 
tions (4) are applied with a very high accuracy. A good 
choice for the phases in Eq. (4) is obtained by es- 
timating the most likely values of these phases which 
will occur during the gate operation (1). We denote this 
choice in the following way 



Ooo 
Oio 



©00 ; 

Oio. 



Ooi 
On 



©01: 

On, 



(8) 



Then, the actual gate is not a perfect gate but it is rather 
the transformation 



(9) 



|00) 


G ^ 


giSoo |oo) ^ 


g»(0oo- 


-Soo) |oo) , 


|01) 




gieoi |oi) — , 


g»(0oi- 


-Soi) |oi) ^ 


|10) 




e*®i°|10) — > 


gj(©10- 


-eio)|io), 


111) 






gj^ gj(0ii- 


-®")|11), 



where we denoted 

i? = eii 



e 



10 



Ooi + e, 



00 , 



(10) 



and we will require i? = tt. The precision of the gate 
is limited by the degree to which G is not perfectly 
predictable from one realisation of the gate to another, 
owing to technical noise and the finite temperature of 
the trapped ions. 



We assume a general initial state 

1 

|*o) = Ca/3|a/3), 

a,/3=0 



(14) 



the minimisation in Eq. (11) runs over all possible initial 

states I'^o) and (•) denotes averaging over other degrees 
of freedom of the system. The imperfection in G can be 
expressed by an error operator E defined through 



G = EGp^rf = E{S^P^) , 



(15) 



where we used the relation = SGpe^f. It follows from 
Eq. (15) that 

E = diag{e'''®'' 



Gperf = diag{e'®™,e' 



^g«eoi gM0io^g«eii|^ (16) 



J&oi gj0io g»®ii 



(17) 



written in the basis {|00), |01), |10), |11)} and we intro- 
duced the notation 

sOc^f} = - . (18) 

Bringing together Eq. (12), (13) and (15) we obtain 

K*perf|*a.t)|' = | (*oiFj.SG|*o) 

= \{^o\PlS{ES^P^)\^o)f 

= K*o|i^|*o)|', (19) 

where we used the commutation of diagonal operators 
and we assumed d — tt in order to accomplish the desired 
gate. Hence the fidelity is 



mm 



E 

a,/3=0 



(20) 



The phases SOap arc random variables and to obtain 
a single fidelity measure it is necessary to average over 
characteristic probability distributions (associated with 
other degrees of freedom of the system). 



C. Fidelity of the phase gate 

In order to calculate the fidelity of the phase gate we 

have to compare an output state of a perfect gate and 
an actual gate. We will define the fidelity as follows 

.F= (minK^perfl^-acOn, (11) 

where the perfect state is produced by the perfect phase 
gate (2) for •& = it, giving 

|*pcrf) =^.|*0), (12) 

while the actual state follows from the actual gate given 
by the transformation (9), that is 

|*act) = SG\^o) . (13) 



D. TT-pulse method 

The general concept in Eq. (l)-(7) of factoring a gen- 
eral phase rotation into single-qubit and two-qubit terms 
can be applied to the error operator E, giving 

E = ZiiSQoo - SOio) Z2i6eoo - SOoi) Psi, , (21) 

where we refer to the notation in Eq. (18), we de- 
noted 61} = 1? — 1? and we dropped a global phase 

(— (58oi — (S8io)/2. We now consider a method which can 
be used to suppress the single-qubit rotations in the error 
operator E, leaving only the part which depends on S^}. 
This method uses a pair of tt pulses on qubits 1 and 2 of 
the form 

R = Ri®R2, (22) 



4 



where 



Ri = |0)i(l|-|l)i(0|, 
R2 = 10)2(11 -11)2(01 . 



(23a) 
(23b) 



We replace the gate sequence SG in Eq. (9) by S'{RG)'^, 
to obtain 



|00) 
101) 
|10) 

111) 



S'(RGf 



gi(5eoo+5eii) |QQ^ ^ 

gi(5eoi+i5eio) |QJ^^ ^ 

gj(5eio+<5eoi) |2^Q^ ^ 
12,9 gi(A-eii+^eoo) 1 ^ 



(24) 



where we used the notation (18), R is defined in Eq. (22), 
G is given by Eq. (1) and 



with 



(25) 



(26) 



where i = 1,2 and a global phase 601 + ©lo has been 
dropped in Eq. (24). Comparing Eq. (9) and (24), it can 
be seen that the joint logical state |11) is rotated during 
the new gate sequence by 2i?, so now we require 2i? = tt. 
This means that in the new sequence (24) each opera- 
tion G will typically be applied for a time half as long as 
in the sequence (9). 



E. Fidelity of the phase gate using the 7r-pulse 
method 



When we use the gate sequence with the w pulses in 
Eq. (24), the actual state of the qubits is 



|*act) = 5' (i?G)2 I fo) 



(27) 



while the perfect state is given again by Eq. (12). Then, 
the fidelity of the gate with the perfect n pulses turns 
out to be 



J^' = ( min 



E 

a, 0=0 



(28) 



where a' = 1 — a and /?' = 1 — fi. We prove this in 
the following (note that Eq. (28) and (37) are identi- 
cal). Whenever SQap ^ 1 and 5<da(3 + SQa'P' ^ SQafj, 
we achieve a significant improvement in the fidelity of 
the phase gate using the tt pulses, that is 1 — !F' <^ 1 — !F. 

The way the 7r-pulsc method works can be brought out 
by writing the gate sequence (24) using the error operator 
expressed in Eq. (21). That is 

S'iRGf = S'iREGp,,if 

= S'[RiR2Zi{^)Z2{^')Ps^Gp,,if, (29) 



where the phases (/) and (f>' are given in Eq. (21), but now 
we do not need to know what they are since we are about 
to prove that they cancel. Expressing the rotations using 
Pauli matrices, we have 

Zj{(f>) = cos{ct>/2)tj+sm{ct>/2){a,)j, (30) 
Rj = iiay)j (31) 



for j 



1,2. Since ay commutes with 1 and anticom- 



mutes with a^, we can write 



RjZj{cj>) = Zj{-cP)Rj. 



(32) 



The operators Zj {((>), Ps^ and Gperf all commute, since 
they arc diagonal, so the rotations Zj{(p) cancel out in 
Eq. (29) because they appear twice in the sequence with 
opposite sign (after using Eq. (32)). Then we get 



S' [RiR2Zii,j>)Z2{<l>')P5^Gp^rif = S' {RPs^Gperif 



The fidelity of the phase gate using the 7r-pulse method 
is 



J^' = ( min 

\l*0> 



(*o| \{GlriR^)\s')n \S'{RG)'] 1*0) 



Using Eq. (29) and (33) this becomes 



mm 

l*o> 



= ( mm 

J*o> 



(34) 



(*0|(Gt^rfi?^) (i?P6,?Gperf)>o) 
(*0|4erf(^^«)'Gperf|*0)f), (35) 

where we applied R^ — R. Let us introduce 

E' = {RPs^f = diag{e*^'', 1, 1, e'^^} . (36) 

Hence, the error operator E' is diagonal so it commutes 
with Gperf, giving 



JT' = (min|(*o|i^'|*o)r 



(37) 



Clearly, when the uncertainties in the phases are such 
that the combination has an uncertainty (5t? small com- 
pared to the imcertainties in the individual phases SOaf} 
in Eq. (16), then the 7r-pulse method offers a significant 
improvement in the fidelity, that is 1 — .F' <C 1 — J-". 

The fidelity of diagonal unitary matrices is discussed 
in Appendix A. We will be interested in the case where 
SQoi ~ — (50io S> (5000,(5011. lu this case the fidelity of 
the phase gate without the n pulses is given by Eq. (19) 
and (A13) as 



.F ~ cos^ [((5001 - (50io)/2] , 



(38) 
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while for the phase gate using the 7r-method Eq. (37) 
and (AlO) give 

T' ~cos^{5d/2). (39) 

It should be recalled that each operation G is typically 
applied for half the time in the gate sequence (24) com- 
pared to the gate sequence (9). Therefore, the quantities 
5Qap in Eq. (38) are not the same as the quantities SQap 
in Eq. (39) as the former are typically twice the size of 
the latter. 

The TT-pulse method ("spin-echo") to suppress an un- 
wanted term in a propagator (evolution operator) is one 
of the standard tools of NMR spectroscopy, which may 
be summarised as the observation that 

[a^Z{ct>)? = a,Z{<p)Z{-(j>)a, = 1 . (40) 

It relies on the assumption that the rotations Z in the er- 
ror operator E in Eq. (21) are the same in two successive 
realisations of the gate operation G in the sequence (24). 
In the terminology of error correction, it takes advantage 

of a known correlation between two errors, thus giving 
an example of the fact that correlation in noise can be 
advantageous. 

F. Imperfect tt pulses 

The TT-pulse method can be a very helpful tool and 
it can strongly suppress the infidelity of a computational 
process. However, it assumes a cancellation of a large 
correlated error, which will only be exact when the cor- 
related error terms arc identical and the tt pulses are 
perfect. Such a cancellation will in practice be sen- 
sitive to imprecisions in the tt pulses and changes in 
the phases 0^/3 . We will next quantify this sensitivity. 

The effect of a change in the phases 6a/3 between 
the two operations G in Eq. (24) can be absorbed into 
the imprecision of the tt pulses. The imprecision in 
the TT pulses can be modelled as a combination of error op- 
erators in the Pauli group. For instance, the imprecision 
in the duration of a tt pulse will give an "over-rotation" 
error of the form 

M - rcos(e/2) -sin(e/2) 
~ [sin(£/2) cos(£/2) 

= cos(e/2)l-isin(£/2)cry, (41) 

where e is the imprecision in the angle through which 
the qubit is rotated in its Bloch sphere, U is the unity 
operator and (jy is the Pauli operator. More generally, 
the combination of the Pauli operators in the dynamics 
will depend on the character of noise in the system and 
may consist of (i) unitary rotations (which may or may 
not be factorisable into a pair of single-qubit rotations) 
or (ii) a non-unitary relaxation. All these cases can be 
treated by using appropriate superpositions or mixtures 
of the Pauli error operators. We analyse next examples 
of noise of both types. 



1. Unitary rotations 

It is convenient to consider a general unitary error 

in the tt pulses as a combination of "over-rotations" in 
Eq. (41) and phase errors corresponding to rotations 
about the z axis of the Bloch sphere (see Eq. (5)). 

First, we will consider pure phase errors. Let us sup- 
pose the first TT pulse applied to ion j adds an additional 
rotation about the z axis by ij, and the second tt pulse 
a rotation hy pjij, where — 1 < < 1. To analyse the ef- 
fects, we replace the two occurrences of R in Eq. (29) by 
RiR2Zi{£i)Z2{e2) and -Rii?2-Z'i(piei)^2(p2e2), respec- 
tively. Using the commutation and anticommutation 
properties as in Eq. (30) (37), the fidelity of the phase 
gate with the imperfect tt pulses is found to be 

= ^min|(*o|^i[(l-pi)£i]^2[(l-P2)£2]i^'|*o)r) 

~ (l-o{j6§ + {l-p)eY)). (42) 

Next, consider pure "over-rotation" errors. To simplify 
the algebra we will assume both qubits have the same 
"over-rotation" angle, but the two different tt pulses can 
have different "over-rotation" angles e and pe. Nothing 
essential is omitted by this simplification. It is found that 

^ ~ [\ - O {e^)]T' , (43) 

in which the coefficient £^ is of order 1 + p^ and J^' is 
the fidelity of the phase gate with perfect tt pulses (see 
Appendix B). 

2. Non-unitary relaxation 

The operator R in Eq. (22) represents perfect tt pulses. 
Let us represent the imperfect tt pulses with the following 
transformation 

4 

Tl{p) = Y,M^pMl (44) 

where p is the density operator describing the state of 
the qubits and are Lindblad operators of the form 

Ml = \/(l-Ci)(l-C2) i?i®i?2, (45a) 

M2 = V(l-Ci)C2 l2, (45b) 

M3 - VCi(l-C2) liOi?2, (45c) 

M4 = VC^liOl2, (45d) 
which satisfy the condition 
4 

Y,MlM^ = t. (46) 

The quantities Qi and C2 are bit-flip error probabilities 
for the qubit 1 and 2. In other words, the operator Mi 
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represents the situation when both qubits are successfully 
Hipped, the operator M2 (M3) corresponds to the failure 
to flip the qubit 2 (qubit 1) and finally the operator M4 
stands for the failure to flip both qubits. 

When using the non-unitary relaxation (44) to model 
errors due to the imperfect tt pulses we describe a state of 
the qubits with a density operator. Then, the complete 
gate sequence reads 

Po GpoG^ 

G [n {SGpaG'^S^)] 
SG [R.{SGpoG^S^)] Gt^t 
7^ (sG [7^ (S'GpoGtS't)] G^S-t 



n 



(47) 



n 



where po = 1^*0) ('I'o I denotes an initial state of 
the qubits. The single-qubit rotations S are given by 
Eq. (3)~(4) with the replacement (8) and we introduced 
additional single-qubit rotations in the form 



S^Si®S2 



with 



Si = 
^2 = 



^1(611 
^2(611 



eoi) 
eio) 



(48) 

(49a) 
(49b) 



where we dropped a global phase 6oi -I- 6io in Eq. (47). 
The state of the qubits after the gate sequence (47) is 
derived in Appendix C. 

The fldelity of the phase gate for an actual state de- 
scribed by a density operator is defined as 



( min (*pcrf |Pact|*porf) 
\ l*o> 



(50) 



where |^porf) denotes the perfect state (12), while pact is 
the actual state (C3) after the gate sequence. The min- 
imisation runs over all possible initial states Fol- 
lowing Appendix C we obtain (C ^ 1) 



(51) 



where T is the fidelity of the phase gate with imperfect 
TT pulses and T' is the fidelity for perfect tt pulses calcu- 
lated in Eq. (28). 

Note that the results (42), (43) and (51) all have a sim- 
ilar basic form, namely that the fidelity T is the product 
of the fidelity T' for perfect tt pulses and a quantity which 
is approximately the fidelity of the tt pulses themselves 
(C ~ £^). The essential point is that the system is not 
especially sensitive to imperfection in the tt pulses, even 
though they are being used to cancel a large correlated 
error. 



III. IONS IN THE MICROTRAPS 

In this section we shall introduce the proposal of Cirac 
and Zoller for the phase gate by pushing trapped ions 
in the microtraps [7, 8]. The basic principle is as fol- 
lows. We consider two ions in two distinct harmonic po- 
tentials (so called microtraps) separated by a distance 
d ~ 1 — 500/Ltm. An internal-state-selective and time- 
dependent force is applied on both ions, such that each 
ion experiences the force when it is in an internal (logical) 
state |1) and no force when it is in an internal (logical) 
state |0) (FIG. 1). 



(a) 



(b) 



(c) 



(d) 



I0> I0> 



I0> 



I1> 



FIG. 1; The picture depicts trajectories of two ions corre- 
sponding to internal states |00), |01), |10), where a time- 
dependent force applies on the ion only when it is in the in- 
ternal state 

We assume that the pushing force is applied on a time 
scale much longer than an oscillation period of the ions in 
the microtraps (adiabatic approximation). We also con- 
sider the displacement of the ions x due to the force to 
be small compared to the separation d between the mi- 
crotraps {x <^ d). Then the Coulomb repulsion energy is 
given by 



id 



1 



1 



1 



d' d-\- x' d- 



X 



(52) 



£ ( X / x\'^ X 

dV^'-d + U '^ + d 



X , 

-d^ ' 



for the internal states {|00), |01), |10), |11)}, where £ = 
g^/47reo. In this series expansion, the constant term pro- 
duces a global phase. The term linear in x produces 
single-qubit rotations, while the term quadratic in x pro- 
duces a two-qubit controUed-phase gate. If the force acts 
for time r, then a phase gate with a phase 1? expressed 
by Eq. (7) is obtained, where 



1 



1 



1 



'0 

\d) dh' 



(53) 



The single-qubit rotations which are produced at 
the same time are rotations through some angles (f>, where 



(54) 
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More precise values of and (j> will be obtained below. 

The advantages of this gate are chiefly: (i) the gate 
time can be short, and (ii) it is found that when a finite 
temperature of the ions is taken into account, the depen- 
dence of the phases on the motional state of the ions is 
small when x <ti d. However the condition x ^ d implies 
that the linear term in the Taylor expansion is large com- 
pared to the quadratic term, and therefore when the gate 
time T is large enough to produce the desired phase = tt 
or 7r/2, the Coulomb repulsion produces a contribution 
to (j) which is large compared to tt (typically (j) ~ SOOtt). 
This means that the singie-qubit phases are highly sen- 
sitive to small changes in the pushing force. Overall, 
a relative insensitivity to motional state is obtained at 
the price of sensitivity to fluctuations in the pushing 
force. In Sec. VI we discuss further methods (in addi- 
tion to the TT-pulse method) to reduce the influence of 
these fluctuations. 



where a = 0,1 refers to the internal ionic state |0) and 
|1), f{t) is a dimensionless time profile and we define 



(58) 



The quantity a is equal (up to the factor l/\/2) to 
the ground state width of a quantum harmonic oscilla- 
tor. In this and subsequent sections it is assumed that 
the force depends on time but not on position. The effect 
of a dependence on position is then discussed in Sec. VII. 



harmonic 
potential 




Coulomb 
interaction 



F^O 




F = 



A. Hamiltonian 

The system of two ions trapped in two separate har- 
monic potentials (microtraps) with a time-dependent and 
internal-state-selective force applied on the ions is de- 
scribed in the semiclassical approach by the Hamiltonian 



Hit) 



a, 13=0 



{a\ (g> 



(55) 



with 



Hap{t) 



2m 2m 

+ ^mLu'^{xa + d/2Y + ^moj'^ix'ij ~ d/2f 

-f (s - - d/2)Fa{t) + (s' -x'p + d/2)Fp{t) 
£ 



(56) 



where the two bare microtraps with a trap frequency to 
are separated by a distance d, m is the ion mass, Xa(t) 
and x'f^{t) are coordinates (trajectories) of ions 1 and 2 
corresponding to their internal states \a = 0,1) i and 
1/3 = 0,1)2, Pa{t) and p'p{t) are momenta of the ions. 
We denoted £ — q/Aneo, where q is the ion charge 
and Eq is the permitivity of vacuum. The parame- 
ters s and s' are associated with the potential which 
the ions experience when they are in their equilibrium 
positions (x — ±d/2) in the microtraps. The quan- 
tity Fa(t) denotes an internal-state-selective and time- 
dependent force which displaces the ion only when it is in 
its internal (logical) state |1). We shall consider the force 
in the form 



Fait) = ahwfit)/a, 



(57) 



FIG. 2: Two ions (representing two qubits) are trapped in 
two distinct harmonic potentials separated by a distance d. 
A state-selective force displaces the trapping potential of 
the ion only when it is in the logical state 1 1) . 



The coordinates Xa and x'^ describe positions of 
the ions with respect to the centre-of-mass of the sys- 
tem. It is convenient to apply a transformation 



Xa - d/2 , 
x' + d/2 , 



(59a) 
(59b) 



where new coordinates refer rather to the equilibrium 
positions of the ions in the microtraps. When we use 
the transformation (59), the Hamiltonian (56) can be 
rewritten in the form 



Hap (t) 



2m 
1 

+ 2' 
1 



+ 







2m 




[{Xa 


•^Oi ) 


[{X'p 





s 



ix'p)^ + 2x'ps'] 



\d + 

where we introduced 



Xait) 

x'pit) 



Fa{t)/muj^ ^aafit), 
F(}it)/mLO^ ^Pafit). 



(60) 



(61a) 
(61b) 



In Eq. (60) the action of the force can be interpreted 
as a displacement of the potential in which each ion is 
trapped and this displacement takes place only when 
the ion is in its internal state |1). The distance scale 
of the displacements is given by the quantity a defined in 
Eq. (58). 
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B. Dynamics 

The dynamics of the ion system is governed by the evo- 
lution operator 



U = Dexp 



h 



Hit') dt' 



to 



(62) 



where the Hamiltonian is given by Eq. (55) and D is 
the Dyson time-ordering operator. It follows from 
Eq. (55) and (62) that 



U\a/3) 



lap), 



where 



e 



a/3 



^J\a0{t')dt'. 



(63) 



(64) 



Thus, the dynamics of the ion system described by 
the Hamiltonian (56) corresponds to a rotation of 
the joint internal state of the ions \a(3) by dynamical 
phases Oq/j. In Sec. IV we shall analyse the phases Qa/3 
in order to apply them in the general discussion of Sec. II. 
However, first we have to understand the motion of 
the ions in the microtraps. 



This means that the time scale t on which the force ap- 
plies has to be much longer than the oscillation period 
Tosc = 2Tr /uj of the ion. The duration of the force r and 
the dimensionless amplitude ^ in Eq. (66) are quantities 
which characterise the force. They have to be chosen such 
that the adiabatic condition (67) and the phase condition 
(i? = TT or 2i? = tt) are satisfied. In the adiabatic approx- 
imation the motion of the ion which is in the internal 
state la) simplifies to the form 



(t) ^ Xa{t) + A{t) , 



(68) 



where the displacement Xa{t) is defined in Eq. (61) and 
the oscillations A{t) arc given by Eq. (D9). 

Now let us consider the case of two ions in two sep- 
arate identical traps. The Coulomb repulsion between 
the ions affects their separation. When two bare micro- 
traps (i.e. microtraps without the ions) are separated by 
a distance d, then the equilibrium distance between two 
ions loaded into the microtraps is bigger by the amount 



Ad= — smh^ <^ - In 
3 16 



where ?7 = 27e/4 and 



r] + l + y^rjirj+i) 



(69) 



C. Motion of the ions 

Once we understand the motion of a single ion in a mi- 
crotrap with an external pushing force, it is relatively 
simple to generalise its behaviour to two ions because 
the presence of a second ion in a separate microtrap 
causes via the Coulomb interaction only small pertur- 
bations (under a certain approximation) to the motion 
of the first ion and vice versa. 

The motion of a single trapped ion is discussed in Ap- 
pendix D and it appears that it consists of three different 
motions. We see from Eq. (D7) that the sloshing mo- 
tion is a fast transitory motion depending on the turn-on 
rate of the force (x cx F). In other words, when we 
turn on the force sufficiently slowly compared to oscil- 
lations of the ion, the contribution of the sloshing part 
to the overall motion will be negligible. Then, we re- 
fer to so called adiabatic approximation. When we as- 
sume the displacement of the potential and oscillations 
of the ion to be of the same order (.t ~ A), then the adia- 
batic approximation can be estimated from the condition 
\6x\ <C min{5, A} as follows 



/* 

J to 



f{t')e''^''dt' 



« f{t) . 



(65) 



where / = df /dt. For a Gaussian time profile of the force 
/(t) = ^e-(*/-)^ (66) 
the adiabatic approximation (65) requires 

g-(a.r/2)2 ^ ^ CJT » 1 . (67) 



Tveomuj'^d^ 



(70) 



gives the ratio between the Coulomb energy and the trap- 
ping energy. For e <C 1 the expression in Eq. (69) sim- 
plifies to Ad « ed/2 but we will be also interested in 
the case when e ^ 1. 

The presence of a second ion also causes a shift in 
the oscillation frequency of the first ion and vice versa. 
In a 3D treatment [8] and taking into account up to 
the quadratic term in the Coulomb energy, the oscilla- 
tion frequency in the longitudinal direction is corrected 
by the factor 



(71) 



with a factor i/l — e/2 for the transversal oscillation fre- 
quency. The higher order Coulomb terms cause the net 
potential at each ion to be anharmonic and we will re- 
tain the third and fourth order in the discussion to follow. 
However, in order to keep the expressions uncluttered, we 
will ignore the distinction between Wosc and wtrap- In any 
case, a fully accurate consideration of this point requires 
a fully quantum treatment [8] of the motion rather than 
the semiclassical one adopted here. 



IV. DYNAMICAL PHASES 

Now we are able to analyse in detail the dynamical 
phases 0„/3 in Eq. (64) which have a key importance 
in the quantum phase gate because they correspond to 
the phases in Eq. (6). We will assume 
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(i) Adiabatic approximation - given by Eq. (67) for 
the Gaussian time profile of the force: ojt 3> 1. 

(ii) Displacements of the potentials and oscillations of 
the ions are small compared to the trap separation: 
ir, A <C (weak force, low temperature). 

(iii) Tight trapping potentials with respect to 
the Coulomb repulsion: e < 1. 

Regarding the third assumption, the analysis simplifies in 
the limit e <C 1, and owing to this some of the analytical 
results to be discussed will only be precise in that limit. 
Namely, e <C 1 allows us to assume the harmonicity of 
the potential what we used in Appendix D. However, 
the gate can still operate with high fidelity when e ^ 1. 

When we now substitute Eq. (60) into Eq. (64) we find 
that we deal with single-particle and two-particle phases. 
Therefore, we write the phases O^/j in the form 



©a/3 = ^Pa -\- + 4>a + (I)'i3 + </'a/3 , 



(72) 



where a, /3 = 0,1 and we introduce so called kinetic 
phases 



= -\j'^^dt', (73b) 



2m 



then the single-particle phases due to the trapping po- 
tential and the force 

1 /•*! 



/ -muP' [{xoc - Xoc)^ - xl, + 2XaS] dt', 

Jto ^ 



(74a) 



\ Imcv' [{x'^ - x'^r - {x'f,r + 2x'0s'] dt', 

(74b) 



and finally the two-particle phases originating from 
the Coulomb interaction between the ions 



1 



^ Jto \d + x'g- X 



■dt'. 



(75) 



We will assume that to < 0, t > and \t\, \to\ r, where 
T is the gate time. 



A. Kinetic phases 

The kinetic phases in Eq. (73) arc associated with 
the kinetic energy of the ions. We can calculate them us- 
ing the expression for the motion of the ions in Eq. (68). 
Taking into account that p = mx and assuming that 
jtoM^I c)o (it means that we turn the force gradually 



on from the zero value and turn it slowly off back to 
the zero) we can write 



+00 

Va = / (xl+2ia^ + A'^^dt 



For the Gaussian time profile (66) we get 



LOT y 8 



and 



II 



COS V' , 



(76) 



(77) 



(78) 



where a = 0, 1 correspond to the internal states |0) and 
^ and r are defined in Eq. (66) and £ with V are 
introduced in Eq. (D9). The contribution (p^ in Eq. (76) 
is a global phase which does not depend on the internal 
state I a) and can be omitted. 

The displacement x scales with the quantity a intro- 
duced in Eq. (58), and by assumption (ii) it is much 
smaller than the ion separation d. Therefore, it is rea- 
sonable to assume that the parameter ^ in Eq. (61) is 
of the order of one or smaller (^ < 1). Then, it follows 
from Eq. (77) and (78) that </?^ < tt, or <C tt when 
the adiabatic condition is very well obeyed. 

The contribution (p^^ vanishes very rapidly with 
the adiabatic approximation (67). Outside the adiabatic 
regime the phase ip^^ would be responsible for a sensitiv- 
ity of the gate to the motional state because scales 
with , where £ is the oscillation energy of the ion. In 
the full quantum approach we would find that the phase 
contribution p^^ corresponds to entanglement between 
motional and internal degrees of freedom of the ions 
which spoils the performance of the gate. 



B. Single-particle potential energy phases 

The phases in Eq. (74) are associated with the trapping 
potential energy and energy introduced to the system via 
the external force F. When we substitute Eq. (68) into 
Eq. (74) we obtain 

^« = \ Z <^^tC^ - VT^aujT^- , (79) 

V o CL 

where a = 0, 1 and we omitted the first term in Eq. (74) 
because it contributes only a global phase. With assump- 
tions (i) and (ii), where a; ~ a, we find (f)a ^ n. We 
obtain the expression for (j>'^ when we replace a with P 
and s with s' in Eq. (79). 
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C. Two-particle interaction phcises 

Finally, we can analyse the two-particle interaction 

phases in Eq. (75). They have a key importance in 
the scheme of the phase gate because they actually pro- 
duce the phase i? in Eq. (6). The phases (pa^ are highly 
nonlinear but we can expand them in a Taylor series us- 
ing assumption (ii). Then we can rewrite Eq. (75) in 
the form 



-A. f 



t oo 

E 

n=0 



d 



dt' . 



(80) 



where 4>af3 arc calculated with precision up to n = 4 in 
Appendix E. The term of order n = contributes only 
with a global phase and can be omitted. The linear term 
(n ~ 1) in Eq. (E6) has a character of a single-qubit phase 
and it is the quadratic term (n = 2) in Eq. (E8) which fi- 
nally produces a two-particle interaction between the ions 
via the Coulomb repulsion. We consider also higher or- 
der terms with n = 3 and n = 4 in Eq. (E9) and (ElO) 
because they contribute with thermal motion contribu- 
tions. In principle it would be sufficient to take into ac- 
count only the cubic term {n = 3) but with the help of 
the TT-pulse method we are able to undo contributions of 
the odd terms. Therefore, we have to consider also the bi- 
quadratic term (n = 4) as a lowest order correction to 
the quadratic term when we apply the 7r-pulsc method. 

Using the results from Appendix E we can calculate 
the phase i? defined by Eq. (7). That is 



1? 



e 



11 



^11 



e 



10 



eol + e, 



00 



(81) 



t>10 - 001 + 000 

2 



(- 

\d 



X 



^ + 6 — 

V2 V^^ 



E2 2 



cos(Ai/)) 



where all single-particle phases lyJa, lyj^, 0^, 
out and only the two-particle phases 0q,/3 
= ipi — il)2 and we introduced the quantity 



(p'f^ cancelled 
contributed, 



(82) 



It follows from assumption (i) that cot ^ 1, from as- 
sumption (iii) that e < 1 and we aim for the phase gate 
with ^ — n OT 2'd = TT. Therefore, we confirm that ^ ^ 1 
as we assumed previously in the discussion after Eq. (78). 

In order to calculate the fidelity of the phase gate 
with ions in the microtraps we need to express the quan- 
tity SOa/3 in Eq. (18). However, first we need to make 
a good estimate for the phases Qa/s in Eq. (8). 

After laser cooling we can describe the oscillation mo- 
tion of the ions with the Boltzmann thermal distribution 
with the probability 



1 



-S/ksT 



(83) 



where fcs is the Boltzmann constant, T is the tempera- 
ture and £ is the oscillation energy of the ion. Further, 
it seems to be realistic to describe the oscillation phases 
^1 and ■02 in Eq. (E4) with a uniform distribution. Then, 
a good estimate for the phases Ga/3 is the average value 



2-K 



d'il;idijj2 



+00 



d£ld£2 (^s,+£2)/kBT 



(84) 



Taking the case that fluctuations in the pushing force 
are negligible, the parameter ^ has a fixed value during 
the gate operation. Then, we can write 



and the gate phase in Eq. (10) is 



-!? = ! 



1 + 



(5) 



2 fe , i2kBT \ 

\f2 hui J 



(85) 



(86) 



where we used Eq. (E12). We require = tt or 2i? = tt to 
realise the phase gate and it follows from assumption (ii) 
that a <^ d (because ^ < 1). Therefore, the phase condi- 
tion reduces to ^ w tt or 20 « tt. It follows from Eq. (82) 
that the phase condition determines directly the gate 
time T (for given values of the other parameters). 

Finally, we can calculate the fidelity of the phase gate 
without the tt pulses given by the sequence (6). When 
we substitute Eq. (Ell) and (E12) into Eq. (20) we get 



?Gr-(5 



(87) 



where the minimisation is carried out in the original 
paper [8] and this result corresponds also to Eq. (38). 
The fidelity of the phase gate with perfect tt pulses de- 
fined by Eq. (28) and corresponding to the sequence (24) 
is 



(88) 



corresponding to Eq. (39). Comparing Eq. (87) and (88) 
it is seen that the 7i"-pulse method improves the fidelity by 
two orders in a/d, and we will see later that a/d 10~^. 



V. PUSHING FORCE 

The state-dependent force was introduced in Eq. (56) 
and Eq. (57), but we have not yet specified how it can 
be realised in practice. In this section we discuss 
the dipole force (originating from a laser beam) as a suit- 
able candidate force. We will consider the dipole force 
in a travelling-wave and a standing-wave configuration. 
The dipole force is always associated with some unwanted 
photon scattering which reduces the fidelity of the gate. 
Quantifying this is the main aim of this section. 
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A. Dipole force 

The dipole force is produced in an intensity gradi- 
ent of the Ught illuminating the atom, typically a laser 
beam which is far detuned from some atomic transition. 
The relevant atomic levels are light-shifted (AC Stark 
shift) creating an additional potential for the particle. 
When the laser is tuned below the atomic frequency (red 
detuning), the particle is driven towards the laser inten- 
sity maximum, i.e. to the minimum of the created poten- 
tial (FIG. 3). For a blue detuning the particle is pushed 
away from the intensity maximum. 



^ Gaussian beam 



e>- 



Dipole 
force , 





|g> 



FIG. 3: Non-resonant laser beam causes AC Stark light 
shifts with a magnitude dependent on the intensity profile 
of the laser beam. The light shifts produce a dipole force 
which pushes the atomic particle towards the minimum of 
the laser-light potential. 

Let us assume a two-level atom (ion) with a lower 
level \g) and an upper level |e) separated by an atomic 
frequency uja- The \g) state corresponds to the logical 
state |1) of the ion while the |e) state is an excited aux- 
iliary state. The logical state |0) is not to be coupled 
to the light (FIG. 4). When we apply a laser beam of 
frequency uj^ the atom-laser Hamiltonian in the dipole 
approximation is 



H 



hujQ 



Ml 



(a+4-(T-) 



(89) 



where ti^ = |e)(e| - \g){gl a+ = |e)(5|, U- = \g){e\, 
n(r,t) is the Rabi frequency on the transition |e) ^ \g) 
and we assumed the laser can be treated as a classical 
light field. In the limit of large detuning A = — ^ 
the standard treatment of this "two- level atom" leads 
to the dipole force (on the atom in the \g) state) 



■ dip 



(90) 



In the travelling-wave configuration we place the ion 
off the axis of the laser beam and use its non-uniform 
intensity profile 



(91) 



where x — xo refers to the position of the ion in the profile 
(cross-section) of the laser beam. We assumed the ex- 
ponential time profile of the laser pulse with a dura- 
tion T and the Gaussian intensity profile with the size 



■|g> 



|0>- 




(a) 



FIG. 4: (a) Logical states |0) and |1) of the ion qubit are rep- 
resented by a ground and excited (metastable) atomic (ionic) 
levels. A non-resonant laser does not couple the state |0). 
(b) Logical states are represented by two sublevels of a ground 
state manifold. Coupling of the |0) state is prevented by 
polarization selection rules, for example using a cr-polarized 
laser. 



of the beam waist w and the radial distance x. We also 
introduced the constant 

67rr/ 12FP 



hck^ hck^w"^ 



(92) 



where F is the line width of the \g) |e) transition, I and 
P are the intensity and power of the laser, c is the speed 
of light and A; = 27r/A is the wave number with the laser 
wavelength A. 





Travelling wave 



Standing wave 



FIG. 5; Scheme of the travelling-wave and the standing-wave 
configuration. The travelling wave is created by a single laser 
beam defiected by mirrors. The standing wave is created by 
two laser beams with different fc-vectors. By adjusting the an- 
gle between the fc-vectors the period of the standing wave 
(node-antinode separation) can be changed. 

In the standing-wave configuration we rather use 
the gradient of the standing field of the laser, giving 



2ino e-(*/-^)%-(-/-)' sm[k{z - zo)] , (93) 



where z — zq corresponds to the position of the ion 
in the standing-wave field. When we now substitute 
Eq. (92) and (93) into Eq. (90) we calculate the pushing 
dipole force on the ions. Comparing the expression of 
the dipole force (90) with the general form of the push- 
ing force in Eq. (57) we can express the dimensionless 
amplitude f , giving 



^trav 
Cstan 



a{x — Xq) 



ujAw 
ak 
'ujA 



|0 |2g-2[(x-xo)/H^ 



(94a) 

rJoP sin[2fc(z - zo)] e-2(-/»)'. (94b) 
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It follows from Eq. (94) that the maximum force is for 
Xq = w/2 in the travelling- wave and for kzQ = n/A in 
the standing-wave configuration. 



B. Photon scattering 

We use a far detuned laser to create a dipole pushing 
force on the ions. Even though the laser is far detuned 
from the \g) |e) transition there is still some photon 
scattering due to a non-zero line width F of the transition. 
A single scattered photon is sufficient to destroy the fi- 
delity of the gate, because scattering only takes place in 
the qubit logical state |1), not in the logical state |0). 
Therefore, a scattering event constitutes a measurement 
of the qubit by the environment. We have to choose pa- 
rameters so that the number of photons scattered during 
the gate is small compared to one (TV <C 1). 

The scattering rate in the two-level atom approxima- 
tion and in the weak-coupling regime (|fi|,r <C |A|) 
is [10] 



R 



r 



2 A2 + (r/2)2 + |f}|2/2 



4 A2 



(95) 



The number of scattered photons during the gate oper- 
ation, i.e. during the period when the laser is on and 
drives the l^) <-> |e) transition, is given by 



N : 



/ R{t')dt', 

J to 



(96) 



where we assume that the laser is gradually switched on 
and off as stated in Eq. (91) and (93), and the com- 
ment below Eq. (75) applies. When we now substitute 
for the coupling constant Q the number of scattered pho- 
tons can be expressed in the travelling- wave configuration 
as 



a, 



and in the standing-wave configuration as 

AT n 2d^^^ 1 



P cos^ {kzo) ' 



where the constants are 



Stan 




(97) 



(98) 



(99a) 



(99b) 



and wc assumed that the distance moved by the ion 
was small compared to the waist size w (period ir/k) 
in the travelling-wave (standing-wave) configuration. 

To calculate the number of scattered photons in 
Eq. (97) and (98) we also needed to express the gate time 



T from Eq. (82), where the phase condition for the gate 
with the TT pulses has to be satisfied. It follows from 
Eq. (86) that 261 w tt, giving 



(100) 



2-K 



2 ' 



and we use the expression in Eq. (94) for the parameter ^. 
Wc also took into account that we apply two successive 
gate pulses. 

When we compare photon scattering for the travelling- 
wave and the standing- wave configuration we get (assum- 
ing Xq ~ UU) 



/W\ 2 



(101) 



where w is the size of the laser beam waist and A is 
the wavelength of the laser light. In practice, we shall 
have w ~ 10~^ m and A ~ 10"'' m, that is A^trav/-^stan ^ 
100. This means that in the travelling-wave configuration 
we have approximately 100-times more photon scatter- 
ing. This is because for a given laser intensity, the smaller 
intensity gradient of the travelling wave (compared to 
the gradient of the standing wave) leads to a smaller 
force, hence a longer gate time and more scattered pho- 
tons. Alternatively, if the force is increased for the trav- 
elling wave case by reducing the detuning, the scattering 
rate goes up as 1/A2 and hence again there are more 
scattered photons. 

The travelling-wave configuration nevertheless has 
an advantage. The ions oscillate in the microtraps due 
to non-zero temperature and experience different instant 
magnitudes of the force in different points within the laser 
beam. In the travelling- wave configuration the excursion 
of the ion thermal motion Aqsc is limited approximately 
by Aqsc ^ w, while in the standing-wave configuration 
it has the stronger constraint Aqsc ^ A (Lamb-Dicke 
regime). Therefore, the latter case requires more laser 
cooling and it is generally more experimentally demand- 
ing. 



C. Contribution of scattering to the infidelity 

We have already discussed in Sec. VB that photon 

scattering appears only when the ion is in the logical state 
|1) and therefore constitutes a measurement of the qubit. 
A measurement will cause a state reduction such that fi- 
delity falls to 1/2 in the worst case. However, depending 
on the atomic structure, photon scattering might also 
be associated with optical pumping of the state from 
the state 1 1) to 1 0) , in which case a single scattered photon 
can reduce the fidelity to zero in the worst case. Making 
the latter (i.e. more cautious) assumption, we therefore 
estimate the influence of photon scattering on the fidelity 
as 



-N 



1-N, 



(102) 
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where TV is the overah number of scattered photons dur- 
ing the phase gate and it only makes sense to consider 
< 1. The number N is given by Eq. (97) and (98) 
depending on the experimental realisation. 



VI. INFLUENCE OF FLUCTUATIONS 
A. Fluctuations in single-qubit phases 

In the original work [8] the single-particle phases were 
regarded as not contributing to loss of fidelity, by the ar- 
gument that as long as their value is known, the correct 
single-qubit rotation operator can be applied to a very 
good precision. However, as we commented in Sec. IV, 
there are present some very large terms in the single- 
qubit phases, so that even a small relative change in 
the value of a parameter may introduce an unknown 
phase change of the order of tt and completely destroy 
the fidelity. For example, when the force is provided by 
a laser beam as discussed in Sec. V, intensity fluctuations 
of the laser cause fluctuations in the pushing force F and 
therefore in its amplitude ^. The single-qubit phases are 
highly sensitive to such changes when a <^ d {or more 
precisely ^ d). 

Bringing together Eq. (72), (77), (79) and (E12) we ob- 
tain the phase required for the single-qubit operation 
in Eq. (4), giving 



$1 



©10 — ©00 

1 



1 s Vs 
+ -- 



d i 1 , ,x 

--+ 0{a/d) 



(103) 



where the first term corresponds to the kinetic phase 
ifl, in Eq. (77) and we dropped (p^^ since it is strongly 
suppressed in the adiabatic limit. The second and 
the third term come from the potential energy phase in 
Eq. (79) and they correspond to the light shift phases 
because the pushing force is represented by the dipole 
force. The next two terms correspond to the linear 
and quadratic contribution from the Taylor expansion 
of the Coulomb energy and 0{a/d) denotes higher order 
terms in this expansion. 

The expression for <I>2 = 8oi — ©oo is the same as $1 
in Eq. (103) except s is replaced by s' and the odd-order 
Coulomb terms in d/a and a/d change sign. Further, 
we will discuss the case when $j can be difi^erent in two 
successive gate operations G in the sequence (24) due to 
the influence of fluctuations of experimental parameters. 

The TT-pulse method goes some way towards alleviat- 
ing the eff'ects of these fluctuations, since then the over- 
all sequence is only sensitive to a change in the single- 
qubit phases between the two successive applications of 
the gate operation G. Let these two successive pushing 
gates be Gi and G2, so the gate sequence for the phase 
gate with the tt pulses is S' {RG2){RGi). Using Eq. (15) 



we define error operators 

where j = 1,2. Then, the gate sequence is 



(104) 



S'{RG2){RGi) ^ y[i?(£;2Gperf)][i?(^lGperf)] . 

(105) 

Without loss of generality, the error operator Ei may be 
taken equal to E in Eq. (21) and we write E2 in the form 



E2 = Zi{A'S>i) Z2{A<^2) Pas E 



(106) 



where A^j are the changes (due to fiuctuations) in $j 
between the two successive applications of the pushing 
gate, that is 



A^j = ($,)^^ - ($,)g^ , 
while Ai? is the change in i9, that is 

A^ = ^G,-^G^- 



(107) 



(108) 



Notice that for A^j = and Ai? = (no fluctuations), 
the sequence (105) corresponds to the sequence (29) and 
both error operators Ei and E2 are equal to E in Eq. (15). 

The imprecision represented by is similar in its 
eff'ects to those of because it reduces the fldelity by an 
amount of the order of (Ad)'^ in the limit Ad <C 1. 

The imprecision represented by A$j can be regarded 
as an imprecision of the second tt pulse and hence is an ex- 
ample of the case which was treated in Sec. II F, namely 
a unitary error in a tt pulse having the form of a rota- 
tion about the z axis of the Bloch sphere. The effect 
on the fldelity is given by Eq. (42) with Sj = A$j and 
Pj = 0. For example, in the case where A^j is caused by 
laser intensity noise, the contribution to the infldelity is 
of order [C{AI/I)]'^ , where A7 is the fluctuation in laser 
intensity between the two applications of the gate opera- 
tion G, and the constant G is of order rf/^a ~ 500. Thus, 
fluctuations of order AI/I ~ 10~^ would be suflicient to 
render the gate useless. 



B. Cancellation of single-qubit phase fluctuations 

We now propose a method to greatly reduce the sen- 
sitivity to fluctuations. We will assume that the fluctua- 
tions A^j are caused by fluctuations in the pushing force, 
and we assume that the force may vary, but s does not. 
This case holds whenever the force comes from a gradient 
of some potential, and the size but not the shape of that 
potential can fluctuate. The analysis therefore applies 
to the case of an optical dipole force with laser intensity 
noise. An analysis with appropriate modifications would 
apply if the source of fluctuations was another parameter 
such as the trapping frequency w hidden in the parameter 
a or the trap distance d. 
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The essence of the idea is to work at a stationary 
value of singlc-qubit phase, i.e. d^^j/d^ = 0. This can 
be done by choosing parameters such that the contri- 
bution from the Ught shift just cancels the contribution 
from the Coulomb energy. The Coulomb contribution 
(the forth term in Eq. (103)) is proportional to the size 
of the force (^ --^ F), while the light shift contribution 
(the third term in Eq. (103)) is proportional to the size 
of light shift itself (,s^ ^ E). Therefore, we have to bal- 
ance the gradient of the light shift with its size. This 
can be done by choosing a convenient distance scale s for 
the laser beams. 

To calculate the optimal value for s, we now have 
to reconsider the analysis, taking force fluctuations into 
account. The effect is that the parameter ^ now de- 
pends on time, where ^ is expressed in Eq. (94) giv- 
ing ^ oc |r2oP oc /. The analysis in previous sections 
is unchanged except that wherever previously we wrote 
^, now we must put a mean value (^), where the details 
of how the averaging takes place will depend on how ^ 
appears in the formulae, and the mean value is taken 
over the gate time r. We will not discuss such details, 
but simply use the formulae, replacing ^" by (^"). To 
estimate the changes in $j due to fluctuations in ^, let 
us calculate the derivative of $j giving 
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(110) 



with a similar expression for d^2/d^ with the replace- 
ments s ^ s' and d ^ —d. 

The pushing force F is derived from a potential which 
we shall call Vp, defined in such a way that for either ion, 
x = (z = 0) is the position occupied by the ion before 
the force is applied. For example, for the travelling-wave 
(TW) and standing-wave (SW) configuration discussed 
in Sec. V we may write 



Y^^-2[{x-xo)/wf^ (TW) 

Vosm^[k{z-zo)], (SW) 



(111) 



where Vq = Vo(t) is a function of time. In what follows 

we use the position variable x but we understand that 
it changes to z in the case of the standing-wave. We in- 
troduced the potential energy into the Hamiltonian (56) 
in the form (s — x ± d/2)F, where after the transforma- 
tion (59) it is (s — x)F. Then, we have 



(s - x)F = Vf{x) = Vf{0) + x^{0) + 0{x'') , 



hence 



s = — 



Vf{0) 



dx 



-riP/AxQ , 
taii{kzo)/2k , 



(TW) 
(SW) 



(112) 



(113) 



where ^^,(0) = ^(0) for TW, ^^(0) - ^ 
SW and we neglected terms of order x"^ (2^) and above. 
The same relations apply for the parameter s' (ion 2). We 
see from Eq. (113) that s and s' can be timed by adjust- 
ing the size of the beam waist w (or standing- wave period 
k) or the position xq {zq) or both. The reference position 
xq (zq) and the waist w can be chosen independently for 
each ion if so desired. The standing- wave period k is tun- 
able without changing the laser wavelength by adjusting 
the angle between the beams forming the standing wave. 

There exists a "sweet spot" or an optimal choice 
for the parameters s and s' which cancels d^i/d^ and 
d^2/d^, namely 

s = -'A + -^[A + eO{a/d)], (114a) 



ed 



s' = +^-^ + ^[A + eO{a/d)], (114b) 

which we obtained from Eq. (109) for d^j/d^ = 0. Typi- 
cally, ed » a so this implies that the parameters s and s' 
have opposite signs. There are two possibilities. Either 
(i) the pushing forces act on the two ions in the same di- 
rection, but the potential energy Vf{0) has opposite sign 
at the two ions, or (ii) the pushing forces act in opposite 
directions and the potential energies have the same sign. 

In the first case (forces in the same direction), 
Eqs. (114) are valid. However, either in the standing- 
wave configuration or in the travelling-wave configura- 
tion (with both laser beams derived from the same laser), 
the light shifts have the same sign at the two ions. There- 
fore, we must have the second case (forces in opposite 
directions) and we have to consider a situation different 
to the one treated up till now. We could eventually con- 
sider the first case if we applied opposite detunings for 
two consequent applications of the gate pulse G. This 
would also directly implement the effect of the 7r-pulse 
method without using any tt pulses [11]. 

With forces in opposite directions, a phase gate is 
still produced with the same pulse timings, but now 
the single-qubit phases are changed. The simplest way to 
analyse this is to make the replacement F/3(t) —Fpit) 
in Eq. (56). It means that we reverse the force on 
ion 2 (compared to ion 1) but we preserve the relation- 
ships given in Eq. (57) and (113). Carrying the analysis 
through, we find that the single-particle potential energy 
phase is 



(115) 



in contrast this with Eq. (79). The two-particle phases 

(j)o.fj arc as described in Appendix E with the replacement 
(3 —j3. Hence we now have 
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so the "sweet spot" is 

s' = +'-^-^[A + eO{a/d)]. (117) 

Note that s and s' still have opposite signs, but now 
(114a) and (117) apply for the situation with light shifts 
at the two ions of the same sign and pushing forces on 
the two ions in opposite directions. Example cases will 
be given in Sec. VIII. 



VII. NON-UNIFORM PUSHING FORCE 

The original proposal [8] assumed that the pushing 
force is uniform. However, we introduced the optical 
dipole force as a suitable candidate for the pushing force, 
where the dipole force originates from a non-uniform in- 
tensity profile of a laser beam. Therefore, the ions expe- 
rience a pushing force which depends on their position. 
In their thermal vibrational motion the ions explore dif- 
ferent positions and as a result the force they experience 
depends on their thermal motion. In this section we will 
analyse this by allowing the dimensionless amplitude of 
the force ^ introduced in Eq. (66) to be a spatial fimc- 
tion rather than a constant (as we have considered up till 
now). 

Using F(x,t) = —dVp{x,t) /dx we obtain from 
Eq. (Ill) the amphtude of the pushing force for 
the travelling-wave configuration 



6rav(2;) = Co 



Xo 



(118) 



where = Ctrav(O) ^^nd xq refers to the ion position in 
cross-section of the laser beam profile. For the standing- 
wave configuration we get 



S,sta.n{z) = C( 



sin[2fc(^o — z)] 
sin(2fc2;o) ' 



(119) 



where = 'Cstan(O) and zq refers to the position of the ion 
in the standing-wave field. 

The non-uniformity of the pushing force introduces 
an additional probabilistic character to the phase &af3 
because the amplitude ^ appears in all expressions for 
single-particle and two-particle phases. Once the value 
of C is not known precisely it introduces a further impre- 
cision to the system and, consequently, a drop of the fi- 
delity. 

The TT-pulse method is again a useful tool because it 
goes some way towards alleviating the effects of the non- 
uniformity in the pushing force. In the full quantum 
approach we would model a thermal state of motion by 
a mixture of vibrational energy eigenstates |n), then each 
state |n) in the mixture acquires a different single-qubit 
phase owing to the different average force it experiences. 
However, as long as the heating rate is low, then the same 
phase will appear in both implementations of the oper- 
ation G and so it is cancelled by the 7r-pulse method. 



Hence, recalling Eq. (36), (37) and (81), we find that 
the fidelity of the phase gate with tt pulses will depend 
only on two-particle phases. To calculate the fidelity we 
need to estimate the difference S-d = between the ac- 
tual and estimated value of the overall phase of the gate. 
In view of the vibrational motion of the ion in the non- 
uniform laser beam profile, a good estimate for the phases 
(compared to Eq. (84)) is now the average value 
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(120) 



where we model the ion's position by a Gaussian proba- 
bility distribution 



P{x) = 



1 



(121) 



The calculation method we have adopted here is slightly 
different from the approach in Section IV. We have not 
analysed the ion trajectories in detail. Owing to the de- 
pendence of the force on position they are now further 
complicated. However, in the limit where the relative 
change in force with position is small, the ion samples 
a small range of force values during its motion. In the adi- 
abatic limit where many vibrational oscillations occur 
during the gate time, we assume the main effect is that 
the size of the force must be replaced by an average value, 
but otherwise the previous calculation of the dynamical 
phases still applies. 

Eq.(121) is the positional distribution for a classical 

particle confined in a harmonic potential at the temper- 
ature T where the standard deviation a is given by 
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(122) 



It is useful to express this in terms of the quantum har- 
monic oscillator energy levels, giving 




(123) 



where (n) « kBT/Tvu) is the mean vibrational number 
and a = ^Jhjmu}. 

Now we have all the tools to calculate the fidelity of 
the phase gate with tt pulses for the non-uniform push- 
ing force. The fidelity depends on the departures of 
the phase from its average (estimated) value, in other 
words on the variance of the distribution of •& values. This 
is calculated in Appendix F, where a complete expression, 
evaluated up to order {a/dY and to [a/wY for the trav- 
elling wave, and to all orders in a/A for the standing- 
wave case, is set out. The most important terms in 
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the travelling-wave case are given by 
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-j Q{2xo/w), (124) 



where Q{y) = 12y^-64y^ +89-34./ y'^+l/y^ and 26»o « tt. 
This expression has been written as a power series ina/w 
and it is sufficiently accurate to give the bclifu iour for all 
the cases considered in Sec. Vlll. The first term corre- 
sponds to the fidelity with the uniform force in Eq. (88). 
However, in practice the other terms almost always dom- 
inate because w < d. The second term vanishes when 
xq = w/2, but away from this position it dominates. 
Therefore, placing the ion at wq = w/2 offers a useful 
improvement in fidelity. 

In the standing-wave case, the expression for the fi- 
delity simplifies considerably when the ion position is 
chosen as kzQ = it /A. This is also the best choice to max- 
imize the fidelity. Restricting to this choice for the ion 
position, the most important terms in the standing-wave 
case are 



32 ^ 



exp 
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(125) 



This reproduces the full expression to good accuracy ex- 
cept when J^' < 0.9. 



VIII. TOTAL INFIDELITY 

We will now bring together all the contributions to 

infidelity which have been discussed, in order to exam- 
ine the performance of the gate in practice. These con- 
tributions are: (i) inaccuracy between actual Oai3 and 
estimated phases Qap, (ii) imprecision of the tt pulses, 
(iii) photon scattering, (iv) laser intensity fluctuations 
and (v) non-uniform character of the pushing force. We 
will use parameters for the calcium ion to give examples. 

Ignoring photon scattering for a moment, it follows 
from Eq. (51) that the fidelity of the phase gate is 



:f={i-ac){i-v').. 



(126) 



where we introduced infidelity V = 1 — !F' and !F' is 
given by Eq. (124) or (125), respectively. The effect of 
laser intensity fluctuations can be absorbed into the term 
C representing imperfection of the tt pulses (see the dis- 
cussion in the last paragraph in Sec. VIA). 

The influence of photon scattering further reduces 
the fidelity. We have not calculated in detail the combi- 
nation of all effects at once, but the expected behaviour 
is that the fldelities multiply, or the infidelities add. In 



the limit of small infidelity these two statements agree. 
Therefore, using Eq. (102) and (126), we estimate the to- 
tal infidelity, including all effects, to be .Ftot = ^'^scat, 
giving 



(127) 



where we assumed all the contributions are small com- 
pared to one (C, V, N ^1) and 26*0 « tt. 



A. Analysis 

We would like to minimize the infidelity while maxi- 
mizing desirable features such as gate speed. The overall 
problem has 11 parameters. Two are associated with 
the design of the ion traps: 

• w - trapping frequency, 

• d - trap separation. 

Four are associated with the laser beam: 

• P - laser power, 

• w {n/k) - size of the waist of the beam (period of 
the standing wave), 

• A - detuning of the laser, 

• xo {zo) - position of the ion in the profile of 
the beam (in the standing- wave field). 

Two describe unavoidable imperfections: 

• T - temperature, 

• C " imperfection of the tt pulses. 

Finally, the last three parameters are associated with 
the ion and the transition which provides the light shift: 

• m - mass of the ion, 

• Aatom - wavelength of a dipole transition, 

• r - linewidth of that transition. 

We will assume the last three parameters are "given", 
i.e. they cannot be varied in seeking the best behaviour. 
We will also assume that C is negligible <C Vtot )- This 
implies that the tt pulses are assumed to be precise, and 
that the laser intensity noise is small. We find that in 
all cases to be considered, the force is only sufficient to 
displace the ions by a small fraction of the trap distance 
{x/d ~ 10~2 - 10~^), and therefore the "sweet spot" 
may be needed to suppress the effect of intensity noise 
sufficiently for this assumption to be valid. 

The position of the ion in the laser beam should be 
chosen sensibly. For the travelling wave, the tempera- 
ture effects (124) are minimized at xq = w/2 (the po- 
sition of maximum force), but photon scattering (97) is 
minimized at xq = w/\f2. Since the fidelity T' depends 
more sharply on xq than the number of scattered pho- 
tons, we will adopt the former choice, i.e. = w/2. For 



the standing wave we chose kzQ = 7r/4 (and xq = 0) for 
the same reason. 

The remaining parameters are lu, d, P, w, k, A and 
T. We will examine a range of values of uj and d, sub- 
ject to technological constraints on what is feasible for 
the trap construction. For the temperature, we keep in 
mind that ion traps in general have three temperature 
regimes which may be relevant: (i) the Doppler limit 
temperature fc^T = hT/2, (ii) cooHng to near the ground 
state (n) « kBT/huj ~ 1, and (iii) ground state cooHng 
(n) <C 1. The regime (i) is achieved by standard Doppler 
cooling techniques, the regime (iii) by the more techni- 
cally demanding sideband cooling, and the intermediate 
case is attractive because it can be achieved by Doppler- 
type cooling on a narrow transition (i.e. using a two- 
photon resonance). Since the gate under consideration 
is relatively insensitive to temperature effects, we will 
avoid assuming ground state cooling, by choosing either 
the Doppler limit temperature or (n) ~ 1. 

Now it only remains to specify the laser parameters P, 
w, k and A. In most conditions higher laser power P acts 
in a beneficial way, so the correct choice is the highest 
power available (exceptions to this rule can occur when 
the laser detuning A becomes of the order of internal 
transition frequencies in the ion). The spot size w (for 
the travelling wave) or the period n/k (for the stand- 
ing wave) has to make a compromise. A smaller value 
increases the force, which reduces the photon scatter- 
ing, but the smaller value also increases the sensitivity 
to temperature through Eq. (124) and (125). 

It is noteworthy that the total fidelity does not depend 
on the laser detuning A (in the weak-coupling regime, 
i.e. \n\ <C |A|). The detuning A can be adjusted to set 
the size of the force (^ oc 1/A) and hence the speed of 
the gate (r oc 1/^^). 

FIG. 6 and FIG. 7 show the total infidelity as a 
function of trapping frequency, for various values of 
trap separation. The laser parameters are set equal to 
the technically undemanding values P = 10 mW and 
w = 4: /im. Results for two cooling regimes are shown: 
either the Doppler limit for calcium on a dipole-allowed 
transition, or cooling to {n) ~ 1. FIG. 8 and FIG. 9 
show results for the more demanding laser parameters 
P — 100 mW and w = 2 /zm. We see immediately that 
good fidelity can be attained, even at the Doppler limit 
temperature. 

Next, let us examine x/d ~ a^/d in order to learn 
the size of the single-qubit phases. Using Eq. (70) and 
(82) we find 




(128) 



where we introduced the characteristic velocity v = 
eV8/{h^) « 1.6(wr)ac, a is the fine structure con- 
stant, c is the speed of light, £ was defined in Eq. (56) 
and 29 w tt. We would like x/d to be small but not too 
small. We therefore choose the gate time r to be small, 
which also gives the fastest possible gate. The limit is set 
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Travelling wave (P=10mW, w=4|im) 




Trap frequency/271 [Hz] 



FIG. 6: The total infidelity given by Eq. (127) as a func- 
tion of the trap frequency for various parameter values for 
the travelling-wave configuration. Results shown are for 
laser parameters P = 10 mW and w — 4 fim. The solid 
lines are for the trap separation d — 100 /im, the dashed 
lines are for d = 10 /im and the dash-dotted lines are for 
d — 1 /im. The thick lines are for the temperature equal 
to the Doppler limit temperature associated with the dipole- 
allowed 451/2 <-* 4Pi/2 transition at A = 397 nm in Ca"*", 
i.e. T = 538 fiK. The thin lines are for fesT — huj, i.e. (n) ~ 1. 
The circles show the "sweet spot" condition (129). 



Standing wave (P=10mW, w=4nm) 
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FIG. 7: Same as FIG. 6, but for the standing-wave configu- 
ration. 



by the adiabatic condition ujt ^ 1. The non-adiabatic 
terms fall off exponentially with (wr)^ so a modest value 
WT = 5 is sufficiently large, hence u = 1.8 x 10^ m • s^^. 
FIG. 6-9 show that the value of usd must be of order 
10^ — 10"^ m • s~^ in order to obtain high fidelity, therefore 
we find a^/d is in the range 2 x 10^^ — 8 x 10^^. This 
implies that in order to keep the single-qubit phase fluc- 
tuations sufficiently small it may be necessary to impose 
the "sweet spot" condition. For a given laser beam size 
and given trap separation, this is a condition on the trap- 
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FIG. 8; Same as FIG. 6, except now the laser parameters are 
P = 100 mW and It) = 2 fim. 



Standing wave (P=100mW, w=2|im) 
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FIG. 9: Same as FIG. 8, but for the standing-wave configu- 
ration. 



ping frequency 



^swcct ~ j\ II 

a\ m\s\ 



(129) 



where we approximated Eq. (114) to the form |s| « ed/4, 
since we find the parameter values are such that the sec- 
ond term in Eq. (114) is a small adjustment. The "sweet 
spot" values of uj are indicated by the circles in FIG. 6-9. 
If the "sweet spot" requirement is imposed, one can no 
longer tune ui so as to minimize the infidelity. 

A further constraint (in addition to the adiabatic ap- 
proximation), which affects the gate rate, is the con- 
dition that the ion is not pushed too far (we will as- 
sume X < w/4 for the travelling wave and x < A/10 for 
the standing wave). Otherwise the ion will move into 
a region where the light shift gradient is substantially re- 
duced, the net result would be a slower gate and increased 



photon scattering. Using Eq. (128) we can write 

COT « I I — , (130) 



Lud 

•^max / 



where aimax = o-S. now denotes the maximum distance 
the ion is allowed to move. We choose Smax = w/A 
(Smax = -^/lO) for the travelling (standing) wave. Tak- 
ing the "sweet spot" value for ujd from Eq. (129), we 
find that Eq. (130) gives the limiting condition on lot 
(i.e. ujT > 5) for the standing-wave case when d > 7 /im 
and for the travelling-wave case when d > 220 /xm at 
w — 2 fim. In general one might expect the larger light 
shift gradient to allow the gate to be faster for the stand- 
ing wave compared to the travelling wave. Owing to 
the previous findings, however this is no longer true when 
the trap separation d is large, i.e. of the order of 50 /im 
or above. 

We find the laser detuning A becomes comparable to 
the fine structure splitting (4Pi/2 ^ 4P3/2) of Ca"^ in 
the case of the travelling wave for P — 100 mW and 
w — 2 /im, but in the other cases the detuning is not 
too large. This would reduce the ratio of pushing force 
size to photon scattering rate, if the qubit is stored in 
the ground state manifold, and therefore result in more 
photon scattering. 

We also find that e ^ 1 for the cases in FIG. 6-9, 
where the fidelity is high. 



B. Discussion 

Although our study has shown that the pushing gate 
in practice does not offer as high a fidelity as was pre- 
dicted in the initial papers which described it, neverthe- 
less comparing this gate with others (which have been 
proposed for ion trap quantum computers), it retains 
a very good combination of allowed temperature and 
gate rate. The 1995 proposal of Cirac and ZoUer [1] 
required ground state cooling (n) <C 1 and required 
the gate rate R to be small compared to rjcu, where 77 
is the Lamb-Dicke parameter [12]. The work of M0lmer 
and S0rensen [13-15] and others has revealed methods 
which allow the requirement on cooling to be relaxed 
to (n) <C l/?7^. The "pushing gate", by contrast, al- 
lows a fast gate rate R — 1/2t ~ w/lO, with simulta- 
neously high temperature. The standing-wave configura- 
tion places the stronger constraint on temperature, but 
even for this case the Lamb-Dicke condition (n) <C 1/77^ 
is sufficient. 

The examples with P — lOmW and w — 4/im show 
that the gate could be demonstrated with modest re- 
quirements on the laser beam power and size. A trap 
separation below approximately 100 /im would be needed 
which is of the order of the smallest values which have 
been reported for ion trap experiments up till now. 

It is reasonable to expect that d = 10 /xm may become 
available in the (near) future. This would permit the in- 
fidehty to fall to Vtot ~ 1 x 10"'^ at the Doppler hmit 



19 



temperature, by using the travelling-wave configuration 
at F = 100 mW, w = 2 jim and u;/27r = 1 MHz, giving 
a gate time of 10/u) = 1.6 fxs. The same laser parameters 
with the standing-wave configuration would allow a gate 
time of 0.3 jis at infidelity Vtot ~ 7 x 10~^. 

In the longer term it would be desirable to have 
a system consisting of a large number (> 10'^) of trap 
pairs, with many gate operations performed simultane- 
ously. With trap separations of the order of microme- 
tres and a gate rate of order 10 MHz, the infidelity of 
order 4 x 10~^ would be available if the laser intensity 
noise were small enough. However reducing the total 
laser power may become an important consideration for 
a large processor. 



IX. CONCLUSION 

The main results of this paper are: 

(i) The single-qubit rotations have contributions from 
the Coulomb repulsion, the light (AC Stark) shift, 
the trapping potential, and the kinetic energy. When 
the parameter regime is such that the two-qubit phase 
is insensitive to ion temperature, then the Coulomb con- 
tribution to the single-qubit phases is much larger than 
TT SOOtt) and proportional to the size of the pushing 
force. This is dangerous because it implies that laser in- 
tensity fluctuations just of the order of 10~^ could be 
sufficient to destroy the fidelity. 

(ii) The dipole force was introduced. It is found that 
the number A'' of scattered photons during the gate varies 
as X/P for wavelength A and laser power P for a laser 
beam focused such that the waist size is proportional 
to the wavelength. Surprisingly, photon scattering is 
independent of the detuning A and the atomic transi- 
tion linewidth F in the limit \^\,T <C |A|, where il is 
the Rabi frequency. The infidelity associated with pho- 
ton scattering increases with the trap separation d and 
the trapping frequency lu, while the infidelity associated 
with finite ion temperature falls with d and lo. The best 
compromise is such that if the separation is small, the fre- 
quency should be large. It is found that very good fidelity 
can be obtained for reasonable values of the parameters 
(see FIG. 6-9). 

(iii) A lack of spatial uniformity in the pushing force 
is found to be a significant consideration when the dis- 
tance scale {w or n/k) of the force profile is smaller than 
the separation d of the ions. In this situation, the effect 
of a finite temperature of the ions is primarily to cause 
fluctuations in the size of the force they experience, and 
hence fluctuations in the phase of the gate. Nevertheless, 
the gate performs well for reasonable values of the pa- 
rameters, and cooling to the motional ground state is 
not needed. 

(iv) Using the tt- pulse method ("spin-echo"), we find 
that a given infidelity (^ or £^) in the tt pulses simply 

reduces the overall fidelity by approximately C or e^. 
An equivalent error occurs if a parameter of the push- 



ing gate changes slightly between two successive appli- 
cations. The most likely source of such fluctuations is 
laser intensity noise. It is shown that the single-qubit 
phases can be made first order insensitive to laser inten- 
sity noise, by choosing the distance scale ("sweet spot") 
of the laser intensity gradient so that the AC Stark light 
shift balances the contribution from Coulomb repulsion. 
This greatly reduces the requirement on laser intensity 
stability, indeed in practical terms it may make the dif- 
ference between a feasible gate and an unfeasible one. 

To conclude overall, the analysis of the points sum- 
marized above was necessary to assess the "pushing 
gate" , and to gain insight into its physical mechanism. 
The analysis which has been presented is sufficiently 
complete to describe an experimental realisation. Even 
though it was found that the fidelity of the gate will not 
be as high as was originally predicted, it still performs 
well in that (i) it does not require ground state cooling 
and (ii) with sufficient laser power neither is cooling to 
the Lamb-Dickc limit required, while the gate time can 
be close to the inverse of trap vibrational frequency, mak- 
ing it a fast gate in comparison with others which have 
been proposed for ion traps. 
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APPENDIX A: FIDELITY OF DIAGONAL 
UNITARY MATRICES 

Suppose a physical operation designed to produce some 
logical operation U instead produces G, where G ^ U. 
The fidelity is then defined to be 

jr = niin|(*o|«7tG|*o)P, (Al) 

l*0> 

where 

1 

a,/3=0 

is an arbitrary two-qubit state denoting \a(3) = \a) (g) \f3). 

Next we will calculate JF for some cases of interest when 
the operator WG is diagonal and unitary. Therefore, let 
us consider a general diagonal operator of two qubits 

U^G = diag{e*®™,e*®",e*®^°,e*®"} . (A3) 

Without a loss of generality we may extract ©oo as 
a global phase and relabel Eq. (A3) as follows 

C/i'G = diag{l,e*^«\e*^i°,e*^"} . (A4) 



20 



When we now substitute Eq. (A4) to (Al) we get 



T = min (1 - 4/) 



(A5) 



with 



f = wx sin2(6'oi/2) + wy sin2(6'io/2) + wz sin2(0ii/2) 

+xy sin2 [(^01 - 0io)/2] + xz sin^ [(^qi - ^ii)/2] 

+yz sin2 [(010 - 0ii)/2] , (A6) 

where to make the expressions less cluttered we used 
the notation 

{w,x,y,z} = {|cooP,|coi|^,|cio|^,|cii|^}. (A7) 

The normalisation w + x + y + z = 1 may be used to 
reduce / to a function of three variables / = f{x, y, z) 
with three parameters ^oi; ^lO) ^ii- 

To find the minimum fidelity the function / must 
be maximised. The general case is straightforward to 
calculate but leads to a complicated expression without 
much useful insight. It is more valuable to consider some 
special cases. 

(i) When one phase is negligible compared to the oth- 
ers (for instance <C 6q\,6\o), using the normalisation 
condition we have 

f{x,y) « {l-x~y)[xsm\eoi/2) + ysm\eiom 

+xys\r?[{en-ew)/2]. (A8) 

Since this is now a function of only two variables x and y 
the maximisation is easier to calculate but still does not 
yield a concise answer. 

(ii) For the case 6 = = ^oi > from Eq. (A8) 
we have 



f{x,y) = {l~x-y){x + y) sin\e/2), 



(A9) 



whose maximum value is /max = (1/4) sin^(^/2) when 
X + y = 1/2. Hence the fidelity is 



JT = 1 - sin2(6'/2) = cos2(6'/2) . 



(AlO) 



(iii) For the case 9 = 9iq = -Oqi > On, from Eq. (A8) 
we have 

f{x,y) = {I - X - y){x + y)sin^(^/2) + xysir? 6 . 

(All) 

The maximum value is on the line x = y, where after 
some algebra we get from the previous expression 

f{x) = 2x sin2(6l/2) [l - 2x sin2(6l/2)] . (A12) 

This function has a maximum value /max = 1/4 when 
e > tt/2, yielding zero fidehty {J^ = 0). For 9 < Tr/2 
the maximum of the function / is when x takes its largest 
value subject to the normalisation constraint 2x < 1, 
giving /max = (1/4) sin^ 9. Hence the fidelity is 



(iv) Finally, if all phases apart from one are negligible 
(for instance 9 = 9oi ^ 9io, 9ii), we have 



f{x)=x{l-x)sm\9/2), 
giving the fidelity = cos^{9/2). 



(A14) 



APPENDIX B: ANALYSIS OF PURE 
"OVER-ROTATION" ERRORS OF tt PULSES 



To study the effect of "over-rotation" in the tt pulses, 
we assume that a single-qubit gate intended to be 
a TT pulse on ion j in fact produces the effect RjMj{ej), 
where M is given in Eq. (41). Therefore, the com- 
plete gate sequence replaces the one in Eq. (24), giving 
S'{RMb)G{RMa)G, where 

Ma = Mi(ei)(g)M2(e2), (Bla) 
Mb = Miipei) (gi M2{pe2) . (Bib) 

Recalling Eq. (34), we find the fidelity is 
min|(*o|<9|*o)|^, where 

Q = [(Gj,,rfi?t)2(5')t] [S'{RMb)G{RMa)G] 
= {gI,,R^)GI^,MbG{RMa)G . (B2) 

This operator, and the fidelity associated with it, does 
not readily simplify. Wc studied it by a combination 
of numerical searches for least fidelity, and by algebraic 
treatment of some less general cases. It was found that 
the case where both ions have the same "over-rotation" 
(ei = £2) does not lead to significantly different re- 
sults from the case where the ions have different "over- 
rotation" . In this slightly restricted case, and considering 
for simplicity a perfect phase gate {G = Gperf), we obtain 



Q 



^+2 



p+i p+i 

p-i -(1+p) 

p-i -(1+p) 

-(1+p) -(1+p) 



s 



l+p'^ - 2ip 







X X 

X X 

. -p^ + i{l + 2p) 



-{2p + p^+i) 





T = cos 



(A13) 



where % = 1 -I- + (1 + i)p. The fidelity obtained from 
this is of the form (a + hp-V cp^)e^, where the coefficients 
a, 5, c depend on the initial state |^o) under consider- 
ation. The main properties of this result are that the 
fidelity is found to be proportional to (not e) and 
there is no special case when p = \ ov p = —\. That 
means the effects of two "over-rotations" by the same or 
by opposite amounts do not cancel. 
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APPENDIX C: FIDELITY OF THE PHASE GATE 
WITH IMPERFECT vr PULSES 

The single-qubit rotations acting on both qubits ex- 
pressed by Eq. (3)-(4) can be written (up to global 
phases) as 



a,/3=0 

= |00)(00|e-*®'"' + |01)(01|e-^®<'i 
+ |10)(10|e-'®i° + |ll)(ll|e 



(CI) 



1011 gi!? 



Similarly, we can rewrite the rotations S in Eq. (48)-(49) 
to the form 

S = J2 \aP){aP\e-'^-^ {e'^Y"''"^ (C2) 

a, 0=0 

= |00)(00| e-'®™e'^ + |01)(01| e"'®" 
+|10)(10| e-'®'° + |11)(11| e"'®" , 

where Sap is the Kronecker delta symbol. The two differ- 
ent single-qubit operations S and S are used here, rather 
than a single operation at the end, because it may be ad- 
vantageous to run an experiment this way, since then im- 
perfections in G can be undone immediately, using a feed- 
forward technique on the hardware which implements the 
gates. 

Even though we introduced different bit-flip error 
probabilities and (^2 for the qubit 1 and 2 in Eq. (45), 

they will be similar in practice (^1 — C2)- Here we will as- 
sume Ci = C2 = Ci where the error is very small {(^ <C 1). 
We note that the overall fidelity of the whole sequence is 
similar in the cases (i = (2 and Ci — C2- 

Following the complete gate sequence (47) we can now 
calculate the actual state of the system (representing 
the evolution of the phase gate with imperfect tt pulses) 
using Eq. (CI) and (C2), that is 

Pact (l-4C)|So)(So| 

+CX:is,)(s,| + o(c^), (C3) 

where 



E 


Ca/3 






'0') (^^i^yP+^c 


\a(3) , 


E 

a,f3 


Ca/3 








\aP'), 


E 


Ca/3 






'0') ^gii5^«/5+'5c.' 


\a'P), 


E 


Ca/3 






'fi) (■g"5-j"/3+'5c,'^j 




E 


Ca/3 






0') ^^i&yi^+^<.0' 





All the summations run over a. /? = 0, 1, we used the no- 
tation a' = 1 — a, /3' = 1 — f3 and applied the approxi- 
mation 



(l + O^^l + mC, C<1- 



(C4) 



In what follows we neglect terms of the second and higher 
orders in ( denoted as 0{(^). 

Substituting Eq. (C3) into Eq. (50) the fidelity of 
the phase gate with the imperfect tt pulses is 



T ~ ( min <^(1-4C) 



+c 
+c 
+c 
+c 



(L—M) gi(^®oo+s@ii) _|_ — M*) e^i^^oi+seio) 



{K~iN*)e 



iiSeoo+SBio) 



+ {K*+iN) e 



.{SBoi+seii) 



(L-iM*) e'('5®'"'+'5®oi) -I- {L*+iM) e'('5®io+50ii) 



(C5) 



where {K,L,M,N} = {coo cgi, cqo cJq, Cn Cgi, Cn cJq} 
and a,f3 = 0,1. The first term in Eq. (C5) is equal 
(up to the factor 1 — 4^) to the fidelity of the gate with 
the perfect tt pulses expressed by Eq. (28). The other 
four terms in Eq. (C5) have always a non-negative con- 
tribution. We can find an initial state (i.e. a combination 
of coefficients Cap) which makes the contribution of these 
four terms equal to zero, thus minimising the expression 
in Eq. (C5). Then the fidelity of the phase gate with 
imperfect tt pulses is as given in Eq. (51). 



APPENDIX D: FORCED CLASSICAL 
HARMONIC OSCILLATOR 

We treat the ion system semiclassically, in that we as- 
sume internal ionic states |0) and |1) but we treat the mo- 
tion of the ions as a motion of classical particles. Then, 
the motion of a single classical particle confined in a har- 
monic potential with an external force is described by 
the Lagrangian 



L{t) = -mx^ — 



(Dl) 



where x{€) denotes the trajectory of the particle, F(t) is 

a time-dependent external force and x =- dx/dt. We can 
rewrite the Lagrangian (Dl) in the form 

L{t) = ^'T^^^ ~ ^rruJ^{x — xf' — ^mw^x^ , (D2) 



where 



a,0 



x{t) = F{t)/muj'^ . 



(D3) 
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It follows from Eq. (Dl) and (D2) that the action 
of the force on the particle corresponds to the time- 
dependent displacement (D3) of the harmonic poten- 
tial in which the particle is confined. The dynamics of 
the particle are given by the Lagrange equation 



where 



An) 
^al3 



+00 

/ 



dt. 



(E2) 



dt dx dx 



(D4) 



Substituting the Lagrangian (D2) into Eq. (D4) we get 
the equation of motion 



x{t)+uPx{t) = F{t)/m, 



(D5) 



which can be solved analytically. The solution splits into 
three terms 



x{t) = x{t) - 5(t) + A{t) , 



(D6) 



where x{t) is the time-dependent displacement of the po- 
tential defined in Eq. (D3), 5{t) is so called sloshing mo- 
tion 



m 



+ cos(a;t) 



/ i{t')sm{ujt')dt' 

J to 



x(t') cos{wt') dt' 



(D7) 



and finally A{t) denotes standard oscillations of the par- 
ticle confined in the harmonic potential 

A(t) = {xo — xo) cos{ujt) + {xo/u) sin(cjt) , (D8) 



where the initial conditions are xq = a;(io), xq = xito) 

dt ' 



and ±0 = ^(io)- We can rewrite Eq. (D8) in the form 



and (jy^^^p = —{£/hd){t — to) is a global phase. Using 

Eq. (68) the motion of ion 1 in the internal state la) 
and ion 2 in the state |/3) can be expressed (when e <C 1) 
as follows 



xc^it) = x„(t) + A(i), 
x'^it) = x'^{t) + A'it), 



(E3a) 
(E3b) 



where the time-dependent displacements Xa (t) and x'p (t) 
are defined by Eq. (62) while the oscillations in the mi- 
crotraps (Eq. (D9)) are 



A(i) 
A'{t) 



2Si 



2^2 



cos{u}t + ipi) , 



cos(wt + V'2) • 



(E4a) 
(E4b) 



The parameters 8\ and £2 are oscillations energies of 
the ions which are in principle different because the ions 
are trapped in two separate and independent potentials. 
Using Eq. (E3) and (E4) we can express Eq. (E2) in 
the form 



An) 



M 



/ 



(a - /3)aJ- +{A- A') 



d 



dt . (E5) 



A{t) = y/2£/muj^ cos{L0t + ijj) , 



(D9) 



where £ is the oscillation energy of the particle in the har- 
monic potential and tp is the initial phase which depends 
on initial conditions. 

We can assiime that the ion experiences a harmonic 
potential only in the hmit e <C 1. Otherwise, Eq. (D6) 
does not describe its motion or describes it only approx- 
imately. 



APPENDIX E: TWO-PARTICLE INTERACTION 
PHASES 

It is sufficient to consider in Eq. (80) with respect to 
assumptions (i)-(iii) only the terms up to n = 4 and 
we choose |to|, |^| ^ 00 (except for global phases) for 
the same reason as it is stated in Eq. (76). Then we can 
have (swapping integration and summation) 



When we now neglect in the adiabatic approxima- 
tion (67) all terms where the factor e~^'^'^^ appears, then 
we can express the linear term (n = 1) from Eq. (E5) as 
follows 



-^(a-/3)ea;r^^ + Qi, (E6) 



where a,/3 = 0, 1, e is defined by Eq. (70), and r are 
introduced in Eq. (66), lu is the oscillation frequency with 
respect to Eq. (71), d is the ion separation with respect 
to Eq. (69), a is defined in Eq. (58) and for global phases 
we used the notation 



Qn 



hd 



to 



A- A' 



dt'. 



(E7) 



The quadratic term (n = 2) reads 



Z>a/3 



'An) 
'al3 



(El) 



(E8) 
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the cubic term (n = 3) reads 



48 



d 



appears in the expression for the fidelity of the gate in 
Eq. (20) and (28). We will find out (in the approxima- 
tion when we ignore fluctuations of the pushing force) 
that the fidelity of the gate is determined by the terms 
of the order n = 3 and n = 4 which represent thermal 
corrections to the motion of the ions. 



— + 2— C0S(V'1 - ll}2) 

nw nui Two 



and finally the biquadratic term (n = 4) reads 



+ Q3 
(E9) 



£1 £2 „ V£i£2 , 1 , ^ 
— + T 2— cos(V'i - W2) 



+ Q4, 
(ElO) 



where £1, £2 and tpi, xp2 are associated with oscillations 
of the ions in the microtraps and they are defined in 
Eq. (E4). 

When we substitute Eq. (82) into Eq. (E6)-(E10) we 
can write the interaction phases 0^/3 in Eq. (75) with 
the precision up to n = 4 in a compact form 



APPENDIX F: NON-UNIFORM PUSHING 
FORCE: FIDELITY 

It follows from Eq. (36) and (37) that the fidelity of 
the phase gate with tt pulses is given for <C 1 as 



(Fl) 



where (•) denotes averaging as stated in Eq.(120). Be- 
fore we calculate the fidelity it is convenient to relabel 
the quantity 6, originally defined in Eq. (82), giving 



(F2) 



where ^0 refers to the amplitude of the force at x = 
(z = 0) and ^ = Ctrav(-T) or ^ = ^stan(z), respectively. 
Carrying out the averaging in Eq. (Fl) the fidelity of 
the phase gate with tt pulses for the non-unitary pushing 
force up to {a/d)^ is found to be 



'pap 



-ia-P)0- ^ 



1 



^ (q - i3fe^ 

a^V2 2 



£1 , £2 ^V£i£2 , 1 , X 
— -h T 2— cos(V'i - W2) 



+3 



£1 £2 „ \/£i£2 , , I S 

— + - 2— COS(V'l - 1p2) 

Two rvaj rko 



(Ell) 



The mean (estimated) value of the phases (j)ai} with re- 
spect to Eq. (84) is 



3na ( £. 1 



6kBT \ 



d\^/6 ^^/2 hoj J 



(E12) 



Now we can easily calculate the difference between the ac- 
tual and the estimated value of the phase (j)a/3 which 



12kBT 
(l2kBT\ 



[f^-(a']+-^(|-6-f2|4) 

J_ ( l2kBT \ 



-(5) m [^«'-<«T-^(^j 



}■ 



(F3) 



where 



+00 



dxC{x)P{x). 



(F4) 



To calculate the quantity ^" which appears in the ex- 
pression for the fidelity, we use Eq. (118) for the case of 
the travelling wave. Then, in the limit r <C 1 , where we 
introduce the coefficient r = 2<t/w « 2ayJ {n)/w, we find 



(F5) 



We will consider also the term of the order of because 

it scales with (a/vj)^ and it is comparable with or larger 
than (a/d)^. Finally, we need to evaluate the coefficients 
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^(") and for n = 2, 4, 6, 8. They are 



^(2; 

^(2 

5(6. 
5(8: 



i?2 



2i?^ 



(F6a) 



3 84 94 

15 , 
— + 18i?2 - 39 , 

45 495 2295 o , 



(F6b) 



(F6c) 



28 

— + 32i?2 - 68 , 



(F6d) 



210 1680 



i?4 B? 



3480 - 2432i?2 + 5127?-* 



where we denoted R = 2x()/w. It is of the order of one 
{R ~ 1) because Xq refers to the position of the ion in 
the laser beam and hence we can assume xq ~ w. 

By similar means, we can calculate ^" for the case of 
the standing wave. Substituting Eq. (119) into Eq. (F4) 
we get 



2sm^{2kzo) 



1-e 



-8(kaf 



COS (4/02:0) 



(F7a) 



8sm'^{2kzo) 



3 _ 4e-8(fe'^)' cos(4fc^o) 



(F7b) 



It is useful to write the exponents in Eq. (F7) in the form 

{kaf « 2rj^{n) , (F8) 

where r] — fca/\/2 is so called Lamb-Dicke parameter, 
k = 27r/A is the wavenumber of the laser light and 
all other constants are defined in Eq. (123). Then, we 
distinguish two extreme cases: (i) Lamb-Dicke regime 
ijf'in) <C 1) and (ii) outside the Lamb-Dicke regime 
{rf'{n) > 1). The former case corresponds to the sit- 
uation where the spatial extent of the vibrational mo- 
tion of the ion is mu ch smaller than the wavelength of 
the laser light {a^J (n) <C A). It means that the ion is cold 
enough to be well localized in the standing-wave and not 
to travel across its nodes and antinodes. In other words, 
in the Lamb-Dicke regime the ion does not experience 
the non-uniform profile of the standing wave. Then, we 
have e-^*^"^)' 1, giving Q and the fidehty (F3) 

reduces to Eq. (88) which describes the case of the uni- 
form pushing force. 

Outside the Lamb-Dicke regime we have the situation 
where e~^'"^^ — > and it leaves us only with the first 
terms in Eqs. (F7). In the expression for the fidelity (F3) 
the first term dominates, giving 



10- rc4 



4^2 L*- J 32sin^(2fczo)' 



(F9) 



The condition for the phase gate gives 26o !v tt and we 
choose kzo = 7r/4, then J^' rs 0.92. 



e = 



so 



10 - 156-^^'='')' cos(4A;^o) 



32 sin^(2A:zo) 
^gg-32(fea)^ cos(8/c2o) - e-^2(fea)^ cos(12fczo) 



(F7c) 



128 sm^{2kxQ) 
^28e-32(fc-)' cos(8yfczo) 



35 - 566-*^'='')' cos(4A;^o) 



-72(fc<T) 



-he" 



-128(/ccr)-^ 



cos(16A;2;o) • 



cos (12/22:0) 
(F7d) 
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